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Abstract 

In quantum chromodynamics (QCD), the vacuum structures are taken to be formed by non- 
perturbative interactions between gluons. The contributions of gluonic interactions can be 
parametrized by gluon condensates and especially the dimension-2 condensate (A'jf) can be nsed to 
generalize the 9 vacunm. In the generalized pictnre, each gauge fixed slice forms a curved surface 
around a hadron and other strongly interacting objects, and the propagation of a glnon has to 
be considered on this cnrved surface. The gluon propagator tnrns ont to be massive due to the 
curvature of the gauge fixed surface implementing the possibility that the massless gauge bosons 
might provide the majority of the hadronic mass and even behave as dark matter in cosmological 
scales. 
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I. INTRODUCTION 


The strong interactions are taken to be mediated by gluons and the most peculiar aspect 
of gluonic interactions is the nonperturbative self-interactions. These nonperturbative self¬ 
interactions generate the problem of definition of vacuum state which is the starting point of 
any calculation by using the properties of quantized fields such as propagators. The vacuum 
states for stroM interactions are known to be classihed according to their topological quan¬ 
tum numbers [l| and those classified states turn out to be mixed through instantons j^. The 
condition that the vacuum state has to be stationary with respect to gauge transformations 
leads to the definition of 6 vacuum resulting from the undetermined angle 6 appearing in 
the eigenvalue of gauge transformation. Because the number of possible values of 9 is inh- 
nite, the gluonic vacuum states can exist in inhnitely different states physically independent 
from each other. However, in ordinary picture, infinitely many vacuum states may gener¬ 
ate inhnitely different vacuum expectation values which lead immediately to contradictory 
results with experimental observations. 

The solution to the problem of inhnitely many vacuum states appeared as a by-product 
of the consideration of in-hadron condensates j4|. The vacuum condensates are dehned 
as the nonzero vacuum expectation values of normal ordered products of held operators and 
are usually taken to be constant with respect to space-time variables. A well-known ehect 
of vacuum condensates is the modihcation of quark and gluon propagators resulting in the 
damping in propagation j^. The damping arises due to background helds residing within 
hadronic scales and therefore the vacuum condensates have to be considered only within the 
scale of hadronic size. In this way the idea of in-hadron condensates has been suggested 
and various examples are discussed with the introduction of maximum wavelength of quarks 
and gluons conhned in hadrons. Although the in-hadron condensate idea resolved many 
problems, there existed intrinsic restrictions concerning the boundary values of condensates 
if the condensate values are taken to be constant inside hadrons. The only breakthrough to 
this restriction is to assume the variation of condensate values resulting in the new picture 

n 

of QCD vacuum [7| composed by union of different gauge fixed slices. 

The new picture of QCD vacuum generalizes the picture of 6 vacuum with coordinate 
dependent 9 values [sj. These coordinate dependences can be parametrized by the values 


of dimension-2 condensates 


Qj and for a given value of dimension-2 condensate the set of 
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points constitute a curved surface around a hadron. Each surface represents one gauge 
hxed slice corresponding to one 6 vacuum and the inhnitely many 6 vacua can be drawn 
as different surfaces around a hadron. Of course this picture can be applied to the more 
complicated situations such as multi-quark states or nuclei and even to the extreme case of 
quark-gluon plasma. In fact, the generalization of 9 vacuum with coordinate dependences 
has been suggested as one possible scenario to explain the strong CP problem in hot 
quark-gluon matter. Recent experiments at RHIC and LHC conhrmed the P-odd charge 
separation effects ll| and the generalized 9 vacuum can be taken as the state where the 
gluonic behaviors have to be analyzed. 

Calculational method for the structures of g aug e fixed slices has been developed by con¬ 


sidering topological spaces of gluonic domains 


12| which can be classihed according to the 


numbers of quarks and antiquarks inside each gluonic domain. For a given topological 
space, the changes of gluonic domains are induced by quark pair creation or annihilation 
and the probability amplitude to have a quark pair is taken to be proportional to the value 
of dimension-2 condensate at that point 13|. When these considerations are combined with 
the amplitude dehned by nonlocal condensate between two quarks, we can deduce functional 
form of the variation of dimension-2 condensate. The surface formed by the points with the 
same value of dimension-2 condensate constitutes one gauge hxed slice and we can draw 
inhnitely many slices by changing the value of dimension-2 condensate. Thus the QCD 
vacuum around a hadron or other structures such as multi-quark states and nuclei can be 
pictured by the union of each gauge slice formed by given value of dimension-2 condensate. 

For a given QCD vacuum, we can try to deduce the Feynman rules and the most important 
quantity is the gluon propagator. In order to deduce the form of gluon propagator we have 
to hx the gauge and the quantization procedure has to be carried out on this gauge hxed 
slice. But now the gauge hxed slices turn out to be curved surfaces and the quantization of 
gluons has to be considered on the curved space. The hrst ehects of curved space are the 
changes of derivatives from the ordinary ones into the covariant ones 1^. These changes 
give rise to additional terms related to the curvature of the space. In the simplest case the 
additional term can be interpreted as carrying out the role of mass term in gluon propagator, 
and thus there appears the possibility that the masslessness of gluons as gauge bosons can 
be changed naturally into the massiveness of hadrons or other strongly interacting objects. 
For cosmological scales these considerations may lead to one possible candidate for dark 
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matter 


id- 


in Sec. II, we review the procedure of generalization of 6 vacuum, and in Sec. Ill, we show 
the structures of gauge fixed slices. In Sec. IV, we deduce the form of gluon propagator in 
curved space, and in Sec. V, we suggest gluons as dark matter candidate. The hnal section 
is devoted to summary and discussions. 


II. GENERALIZED 9 VACUUM 


For atomic structure, the composing elements are nucleus and electrons and these ele¬ 
ments have dehnite masses with nonrelativistic movements. The vacuum state for an atom 
is taken to be a state in which nothing exists. From this vacuum state we can count the 
number of electrons and an atom is formed with these dehnite number of electrons and 
nucleus. However, as we approach the nucleus it becomes impossible to count the number 
of electrons because electron-positron pairs are continuously created and annihilated as the 
energy of the probing particle goes beyond 1 MeV. Then the vacuum state for a nucleus 
cannot be a state of nothing. It is a state of somethings that are created and annihilated 
continuously. In order to describe physical quantities systematically we need to arrange 
the creation and annihilation of particles in a dehnite manner and the resulting method is 
known as normal ordering. In quantum held theories, the method of normal ordering is fully 
exploited to dehne physical quantities. 

The situation changes once again when we consider nonperturbative self-interactions of 
gluons. The normal ordering method is useful only when we can count the number of 
particles in a given state. But due to the nonperturbative self-interactions the number of 
gluons in a state cannot be counted appropriately. Therefore vacuum expectation values of 
normal ordered held operators do not vanish for gluonic states and these are parametrized 
as vacuum condensates. The original idea of condensate was applied only to gauge invariant 
forms of held operators and the held operators were classihed according to their dimensions. 
However, after three loop calculations had been carried out [l6|, there appeared the necessity 
of introducing dimension-2 condensate which depended on the choice of gauge. If something 
is dependent on the choice of gauge, it can be used to hx the gauge. Since the quantization 
procedure of gluons has to be carried out with gauge hxing, the slice with given value 
of dimension-2 condensate can be the appropriate place where the quantization of gluons 
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could be performed. The collections of these gauge slices are intimately related to the 
generalization of 6 vacuum. 

The dehnition of 9 vacuum starts from the instanton solution to minimize the Yang-Mills 
action satisfying 

F^u = g[A^^ AJ\ = 0. (1) 

To induce the form of A^ satisfying these equations, we consider the transformation prop¬ 
erties of Afj, at long distances and get 

A^{x) = ( 2 ) 

where g{x) represents the matrix element of the group of gauge transformations. In Eu¬ 
clidean space with imaginary time, the set of points at inhnity in four dimensional space 
forms a large sphere and the above relation generates mappings between the points on 
and the group elements g{x). These mappings classify the gauge helds A^[x) into inhnitely 
many independent components. For a group element g'i(x), the gauge field Ai(x) for a 
vacuum is given by 

Ai(x) = ^h^(x)V^i(x). (3) 

Representing the mapping corresponding to n times repetition of 5 'i(x) as 

^n(x) = [c/l(x)]^, (4) 


the gauge held given by 

An(x) = ^■^(x)V^„(x) (5) 

constitutes another vacuum state locally. When we write the wave functional for this vacuum 
state as t/’n[A], the other vacuum state '^„+i[A] at A„_|_i(x) can exist and these states can 
be mixed through tunnelling. Then the true vacuum state can be written as 

'l'|Al = (6) 

n 

and the coefficients Cn can be determined from the condition that 4/[A] has to be stationary 
with respect to gauge transformations. The hnal result is 

>l<»[A]=^e“V„[A], (7) 

n 

and this state is called the 6 vacuum. 
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Now the parameter 6 remains as a free parameter. For two different 9 and 9' the transition 
amplitnde becomes 


{9'\ 




\9) 


6{9 - 9') Y, e"*"' J [DA^ ■■■], exp{- J 


]^exp{- / ••■)}• 


( 8 ) 


The exponential factor in front of the fnnctional integral represents the change of vacnnm 
winding nnmber and this factor can be inserted into the fnnctional integral by changing the 
Lagrangian density as 


C = tr 




( 9 ) 


in Minkowski space. The second term breaks the P and CP invariances and these effects 
have been checked by the measnrement of nentron electric dipole moment resnlting in the 
limit 9 < 3 X 10“^°. However, in recent experiments at RHIC and LHC, local P-odd effects 
are observed implementing the possibility that topological flnctnations exist in hot qnark- 
glnon matter. Then there exist dehnite differences between the glnonic states of a nentron 
and hot qnark-glnon matter. Becanse the differences are generated by instantons, we need 
to analyze the relationships between the instanton contribntions and the strnctnres of QCD 
vacnnm. These relationships have been checked by considering the correlation fnnction 
on a lattice in one approach, and by nsing the instanton shape recognition procednre on 
the lattice in another approach. From the comparison of two-point correlation fnnction 
with the lattice glnon propagator, 0(p) correction to glnon propagator has been conhrmed 
indicating the existence of dimension-2 condensate 17(]. By using instanton shape 

recognition procedure ISj we can also conhrm the presence of dimension-2 condensate and 
the expected values are well explained in an instanton-gas or an instanton-liquid picture. 

For topologically non-trivial vacuum states, the strong interaction Lagrangian is given by 


r = _ 

^ 4 Ml" “ 


327r2 


f 


( 10 ) 


The ^-term characterizes the vacuum state and the observations of parity-odd domains 
localized in space and time strongly support the spacetime dependence of 0 as 6* = 9(x.,t). 
In case of varying 9 values, the vacuum domains can be classihed by the conditions 0(x, t) = 
9i with hxed values of 9i. The classihed points constitute surfaces and according to the 
movements of quarks the forms of surfaces will change. For example, the points with hxed 
9iS for 6 quark conhguration can be drawn as in Fig. [T] 1^. This conhguration can be 
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made in the scattering processes of proton-proton collisions. Since the different surfaces 
correspond to different values of 6, we have to consider the strong interaction vacuum as 
formed by the union of different slices in case of spacetime dependent 6. For a given slice, 
the instanton contribution can be estimated by the value of dimension-2 condensate 
Then we can assign another value of to another slice. These conhgurations can be 
represented as 

{^l)e = Ce ( 11 ) 

with hxed C 0 for each slice. For different slices we can write 

(Aj) = C(x,t), (12) 

and the spacetime dependence of 9 is transferred into that of (^^). The whole QCD vacuum 
is now composed of these slices establishing the generalized 6 vacuum j^. 


III. STRUCTURE OF GAUGE FIXED SLICES 


In general relativity, gauge invariance stems from the independence of the held equations 
on the choice of coordinate system. The Einstein held equations depend only on metric 
tensor and energy-momentum of the system, and therefore the choice of a coordinate system 
does not ahect the equations of motion. However, in order to solve the equations of motion 
it is quite important to choose the coordinate system and this choice of particular system 
corresponds to the hxing of particular gauge 2l|. For quantum held theories, the main 
variables are changed into the held variables at each space-time point, and since the held 
variables generally include internal degrees of freedom the gauge hxing condition becomes 
functional relation specifying the surface on which each gauge orbit has unique position. 
The gauge orbit is obtained by connecting the held variables that can be related by gauge 
transformations. Thus the determination of functional relation for gauge hxing is a dual 


process with respect to the following up the gauge orbit 


22 |. 


For the space of gluonic held A„, the condition for gauge hxing can be written as 


/(A,.) = 0. 


(13) 


Traditionally the functional form of / had been taken to be linear in A^ because of simplicity 
and the corresponding surface is taken to be a plane in the held space. However, now 
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we have the generalized 6 vacuum dehned by the union of gauge slices given in Eq. flT^ 
where the functional form is quadratic. The quadratic form originates from the existence 
of dimension-2 condensate due to the nonperturbative self-interactions of gluons. These 
gluonic self-interactions induce that the gluonic states cannot be classihed according to the 
number of gluons and the unique possibility to dehne the gluonic state is to consider the 
state of inhnite number of gluons. Let this state be represented by |r2) and then we have 

Al{x)\n) = (5\n) (14) 

because the gluon number is not changed by one creation or by one annihilation. We get 
the covariant relation 

{Sl\Al{x)A^^{x)\n) = a^ (15) 

with inner producting each side. This condition is equivalent to the gauge slice condition 
in Eq. fll2p if we take the state |r2) as one vacuum state. For different a values we have 
independent vacuum states represented by different values of dimension-2 condensate. In 
ordinary flat space, these vacuum states might be overlapped, however, because the vacuum 
condensates have to be considered only around hadronic structure as indicated by in-hadron 
condensates, the QCD vacuum is taken to be the union of curved slices hxed by the value 
of dimension-2 condensate. 

In order to hgure out the structure of gauge hxed slices, we need to classify the gluonic 
domains formed by slices with given values of dimension-2 condensate. Because the gluonic 
domains are made around quarks and antiquarks which behave as sources and sinks of the 
gluons, we can classify the domains by the numbers of quarks and antiquarks in each domain. 
If we represent the gluonic domain around a quarks and b antiquarks as then meson 
domain is represented by i?ij, and baryon domain corresponds to R^.^. In this representation 
we can include arbitrary number of quarks and antiquarks such as i? 2 , 3 ; however, since the 
hnal states of strong interactions are observed as hadrons we need to categorize the gluonic 
domains with dehnite rules of combinations. Two gluonic domains are united when one 
quark-antiquark pair annihilates and conversely one gluonic domain can be divided into two 
domains by one quark-antiquark pair creation. Therefore the union and the intersection 
operations can be dehned into the sets of gluonic domains, which can be formulated as 


follows 23|: 

• Gluonic domains are open sets. 


7 



• The union of gluonic domains is a glunoic domain. 

• The intersection between a connected gluonic domain and disconnected gluonic do¬ 
mains is the reverse operation of the union. 

These are just the conditions needed to construct the topological spaces of gluonic domains. 
The constructed topological space encompassing i baryons and j antibaryons can be repre¬ 
sented as 


r,j = {0, 


(16) 


For the constructed topological spaces, we may introduce a systematic measure to deduce 
the forms of the gauge hxed slices dehned by the values of dimension-2 condensate. The 
measure is related to the nature of QCD vacuum and so it has to be represented by quantities 
characterizing the vacuum. Because the characteristics of vacuum states have been expressed 
as various condensates, it is plausible to compare two fundamental condensates to deduce 
the form of gauge hxed slice. The two fundamental condensates can be chosen as the quark 
condensate and the dimension-2 gluon condensate. The nonlocal quark condensate is dehned 
by |24| 

(: g(x)f/(x,0)g(0) :) = (: g(0)g(0) :)Q{x^), (17) 


where t/(x, 0) represents the connection through the gluonic domain. In the gluonic domain, 
the gluonic ehects are parametrized by the value of and this value changes from point 
to point. The comparison with quark condensate can be made possible by assuming 


(: q{x)U{x,y)Al{y)A^^{y)U{y,Q)q{Q) :) oc (: q{x)U{x,y)q{y)q{y)U{y,Q)q{Q) :), (18) 

where the hnal result can be expressed as products of two Q’s. To derive the functional 
form of Q{x^) we consider a measure Tt{Q). In general, we can assume that the value of Q 
decreases as the distance x between the quarks increases. This assumption can be stated as 


91I((5) increases as Q decreases. (19) 

Another condition can be written down by comparing both sides of Eq. (fT8|) resulting in 


TliQM = m{Qi) + m{Q2). 


( 20 ) 


Then these two conditions lead to the solution 


m{Q) = -k\n^ 

Vo 


( 21 ) 



where k is an appropriate parameter and Qq is a normalization constant. When we try to 
represent the measure Vyt{Q) as a metric function, we may use the distance function 

d(x,y) = |x-y|^ (22) 


with X and y representing the positions of the quark pair. We need to sum over the contri¬ 
butions from different values of u and the full amplitude becomes 

Q = Qoexp ^ . (23) 

The weight factor F{i/) has been introduced to account for possible different contributions 
from different zz’s, and the variable r is 


?" = y |x — y| 


(24) 


with £ being a scale parameter. For the case of equal weight F{v) = 1 and for a = 2, we get 

f 1 — r 1 

Q = Qoexp <; -- ^ \. (25) 


k In r j 

Perturbative local effects can be included by substituting r^Q for Q, and then Q becomes 


Qo f Ir^-r 


Introducing the notation 


Qij 


QiFj) _ Cdr* -r^l) 


Qo 


Qo 


(26) 


(27) 


we can write the value of {Ai) at the point x as 


6 I 6 

i^i )=^0 n 1 n Qjk ■ n 

rj,rk^ri rj,rk ra,r^^rj,rk 


(28) 


2=1 


for the case of a 6-quark domain and several slices are shown in Fig. [T] For a meson domain 
1 , curved slices are shown in Fig. [2l 


IV. GLUON PROPAGATOR ON CURVED GAUGE SLICE 

Before we discuss the gluon propagation in curved space, it is valuable to summarize the 
changes of ordinary gluon propagator when the dimension-2 condensate exists. For 
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simplicity we assume that the space is flat and the value of is taken to be a constant. 
Because the condensate value is taken to be constant, we may decompose the gluon held as 


.4“(,t) = A“+.4“(,t), 


(29) 


where Af, contributes to the condensate as the zero-momentum mode. Then, the effective 

r—I 

Lagrangian density for gluons can be written as 251 


Ca = - + ^AHx)A‘^(x) 


m: 


with 


GUx) = d.Atix) - d.AAx) + gr^AAx)At{x). 


The gluon mass me is defined by the condensate value 


/ I A rt A I \ r 

(vac |A^A„| vac) = -ai^Aab-^, 


(30) 


(31) 


(32) 


and the existence of mass term indicates that the polarization tensor for a gluon has the 
contribution 


frab 

/i.z/(mass 


) “ [s^lu -p-J 


m: 


G- 


(33) 


If we define the scalar function by 


Kbik) = - k^k^) 


(34) 


then n(A:^) can be written as 


ml U2{k^) 

The first term represents the mass effects originating from the vacuum condensate corre¬ 
sponding to the Schwinger mechanism, and the second term appears as the next corrections. 
The form of 112(A:^) is determined by the relation 



U{k^ 



where 


with 



(36) 


(37) 


j;{x) = sr*” [A",(x)d,A‘^(x) - 2At(x)d,Ai(x) - A((i)aaA‘(,T)] 


+ A’’<x)A'Hx)Ai(x) 


(38) 
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By putting the condition in Eq. 


i J d^x ^vac 
and then we get 


T 


KWtM 


and using the relation in Eq. 

k^k’' 


vac > = 


— 


9^ 


k^ 


■^ab 


, we obtain 


■^mQk"^ 


_F-m^-n2(P)_ 


Mk^) = 


^rriQk'^ 


fc2 


m 


2 • 

G 


With these results we can write the form of transverse gluon propagator as 

D^TMik) = ^ 


A;2 




(39) 


(40) 


(41) 


3 k’^—rtiQ 


where we can see that there exists no pole and so the gluons cannot propagate freely reaching 
far away from the source. This damping of gluon propagator originates from the constant 
dimension-2 condensate generating dynamical mass of gluon. Phenomenologically the dy¬ 
namical gluon mass is determined to be about 500 ~ 600 MeV 26|], however, these values 


are too large to be considered as the value of just one gluon because there exist so many 
gluons in a hadron. Moreover, the gluonic structure of hadron is not uniform as assumed in 
Eq. fl3^ . so that we need to quantize the gluonic held on the curved gauge slices as shown 
in Fig. [Hand in Fig. [2J 

The gauge slices are dehned by the Eq. (IT5|l with different a values and they turn out 
to be curved surfaces around quarks and antiquarks forming hadronic structure. In order 
to account for the curvature we need to replace the ordinary derivative by the covariant 
derivative. The covariant derivative is dehned as 

dKy 




■pAC A 


(42) 


where Tjj^ represents the affine connection which relates the ordinary coordinates to the 
freely falling coordinates. When TJj^ = T^;^, the diherence between two covariant derivatives 
is given by 

(43) 


D A -DA 


dx'' 


which implies that covariant curl is the same as ordinary curl and there appears no curvature 
ehect in hrst order derivatives. However, the Lagrangian which is used in the derivation of 
the gluon propagator contains the second order derivatives and the diherence between two 
second order derivatives becomes [27j 




(44) 
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where the Riemann-Christoffel curvature tensor is dehned as 


dx^ 


(lUK. 


OT^ 

I t^T] -per 

dx'' 


pTJ po- 
■*' flK,^ VT}' 


The Ricci tensor is given by 


p _ pt^ 


(46) 


(46) 


and the curvature scalar is dehned as R = 

Now the free Lagrangian for gluon propagator in curved space can be written as 

ic = -i / i'x (D„A„ - Z)„A„). {D‘‘A‘- - D-A"). 


(47) 


In order to carry out the functional integration we have to arrange the differential operators 
into the form 

LD = \j d^x D^) A,, (48) 

where the covariant differential operators do not commute as in Eg. (1441) . To compare with 
ordinary gluon propagator in hat space we rearrange the diherential operators and get 

1 


Lb — 


d^x A, 


(49) 


with additional Ricci tensor R^^. The gluon propagator is obtained as the inverse matrix 
of the tensor that appears between the two gluonic helds. In hat space, the diherential 
operators can be replaced by momentum variables and the inverse matrix can be calculated 
easily. However, in curved space, the appearance of Ricci tensor makes it non-trivial to 
calculate the inverse matrix and we need to introduce additional assumptions to deduce the 
form of gluon propagator. The simplest case is that of hat space in which the Ricci tensor 
vanishes and we can get the ordinary form of gluon propagator. The next case is that of 


maximally symmetric space 


28| where we can put 




(60) 


with appropriate parameter nR. This parameter is dehnitely related to the curvature scalar 
of the space and is obtained by considering free propagation of massless gluon as a gauge 
boson. Then the Lagrangian becomes 


Ru) = 2 / d X A^ 


+ m^)g>^'' - A^, 


(51) 
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and the gluon propagator can be obtained as 


(52) 


> k^-m^ 

which is just the propagator form of a particle with momentum k and mass m. 

In the usual process of dynamical mass generation of gluon propagator, the polarization 
tensors are calculated by estimating various loop contributions with appropriate renormal¬ 
ization constants. Therefore the generated mass always has some uncertainty related to the 
renormalization process and the loop expansions turn out to be meaningless for nonpertur- 
bative gluonic interactions. Moreover the determined mass is too large to be considered as a 
correction to massless gluon. Instead of these difficulties, it is quite natural to get the mass 
effects of gluons from the curvature of gauge hxed slices. 


V. GLUONS AS DARK MATTER 


The curved gauge slices can be drawn not only for hadrons such as mesons and baryons 
but also for any strongly interacting systems of quarks and gluons such as multiquark states 
and nuclei and even for the system of quark-gluon plasma. Because gluons are taken to be 
propagating along these gauge hxed slices, the curvature effects are essential in estimating 
the gluonic contributions to strongly interacting systems. For perturbative approximations, 
gluon propagators are assumed to be in the same form as that of photons except for color 
factors, and therefore the main features in asymptotically free region resembles those of 


Abelian gauge theory 


29| . But in the infrared region, that is, in the conhning region. 


nonperturbative effects play important roles in changing the vacuum states with nonvan¬ 
ishing condensate values. These changed vacuum states can be classihed according to the 
conhguration of the source quarks. 

For stable hadrons, the conhgurations of quarks are very simple. In a meson, only one 
quark-antiquark pair exists and we can easily draw the gauge slices as in Fig. [2l However, 
when we want to probe the internal structure of the meson, the main scattering processes 
are related with the quarks and the asymptotically free gluons near the quarks leaving aside 
the whole bunch of gluons described as gauge slices. If the probing particle has the enough 
energy to approach the asymptotically free region, it can create quark-antiquark pairs and 
the created quarks and antiquarks change the structure of gauge slices and form the hnal 
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states with new gauge slices. Because the hnal states are counted by hadrons instead of 
quarks and gluons, the nonperturbative gluons described as gauge slices behave just like the 
dark matter. They were not probed but contributed to the total mass of the hnal states. 

In a baryon, three quarks have been introduced to account for the quantum numbers but 
their masses could not be determined in a unique way. The simplest way to determine the 
quark mass is to divide the total mass of a baryon by 3 and the determined mass is named as 


constituent quark mass 


30l |. Although the lowest baryon states, the proton and the neutron. 


have nearly the same mass, the other excited baryon states have quite different masses and 
therefore the value of constituent quark mass has to be different from baryon to baryon. 
These difficulties have lead to the idea of dynamical quark mass |3l| including the effects 
of gluonic contributions, however, in this case also, the probed physical mass is assigned to 
the quark leaving the gluons as dark matter. This situation is not changed in scattering 
processes. In scattering processes, baryons are treated as composed of partons which carry 
fractional momentum of the baryon and the hrst approximations were made in such a way 


that the quarks were taken to be the partons without considering the roles o 


approximations resulted in various forms of quark distribution functions 32| but could not 


gluons. These 


correspond to the experimental data. The corrections to the hrst approximations had been 
made in the direction of introducing sea quarks and there appeared various models with 


different forms of valence quark and sea quark distribution functions 


33| . Because the 


probing particles for the internal structures of baryons were leptons such as electrons or 
neutrinos, the main processes occurred through the exchange of photons or weak gauge 
bosons. Since the gluons did not interact directly with these gauge bosons, it was natural 
to introduce sea quarks instead of gluons to account for the hnal states. However, the total 
momentum carried by charged partons turned out to be only one half of the total momentum 
of the probed baryon. The neutral partons had to be introduced playing the role of dark 
matter. 

The hadronic masses are basically measured by the trajectories in magnetic helds and 
then compared with other scattering data. The trajectories in magnetic helds are determined 
by the charge and the momentum of the hadron. The charge of a hadron is taken to be the 
sum of charges carried by valence quarks, and the momentum of the hadron is measured by 
the velocity of the center of mass of the hadron and the total rest mass of the hadron. But 
the total rest mass of a hadron is not given by the sum of rest masses of the valence quarks. 
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and in fact, consists mostly of the intrinsically massless glnons moving aronnd. In order to 
dednce the rest masses of valence qnarks, we need to get rid of strong interactions and by 


considering electromagnetic and wea. 
qnark and 4.8 MeV for down qnark 


i interactions we can get the values of 2.3 MeV for np 


34| . If we snm np these valnes, we have 9.4 MeV for a 


proton and 11.9 MeV for a nentron. These valnes are just 1 ~ 1.3% of the total masses of 
proton and neutron. Then from the viewpoint of measurement more than 98% of nucleon 
masses are dark matter which is dehnitely composed of gluons. 

To establish the standpoint of gluons as dark matter, we need to check whether they can 
have survived from the Big Bang and can exist in galactic space in a nonbaryonic state. It is 
helpful to use the arguments of maximum wavelength of conhned gluons applied to induce 
the idea of in-hadron condensates. It has been suggested that gluons conhned in a stro ngly 


interacting system have a maximum wavelength of the order of the conhnement scale 


^1 


because they are conhned within that system. Thus gluons conhned in a hadron have the 
maximum wavelength 


\ 1 f m ~ 

/'max — -L — 


200 MeV 

and those conhned in a nucleus of diameter of 10 fm have 

Amax ~ 10 fm ~ 


20 MeV 


(53) 


(54) 


During the period of big-bang nucleosynthesis only light nuclei such as D, ^He, ^He, and 
^Li have been formed, and so the maximum wavelength of the gluons conhned in such light 
elements can be estimated as 


Amnx ~ 2 fm ~ 


100 MeV 


(55) 


This implies that of the gluons created at the Big Bang only those with energies greater 
than 100 MeV could be captured inside baryons or light nuclei. Therefore the gluons with 
smaller energies had to be left as sterile gluons if there existed no larger structure formed 
by strong interactions. The heavier nuclei had been formed later in stars and in supernova 
explosions and the gluons with energies between 20 MeV and 100 MeV could be captured 
in those elements. Then how about the gluons with energies less than 20 MeV? Are there 
any larger structures than the heavy nuclei interacting strongly through exchanges of low 
energy gluons? 
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One possible candidate to accommodate the long wavelength gluons is the neutron star. 
Neutron stars are known to be formed by gravitational collapse after supernova explosions 
and the average density is approximately equal to the density of nucleons. This means that 
the neutrons inside neutron star cannot exist as independent neutrons but have to be in 
the form of overlapped conhgurations. Thus the outer part of a neutron star is taken to 
be in the state of solids or liquids composed of neutrons and the inner part in the state of 


quark-gluon plasma 


In the quark-gluon plasma region, the average distance between 


quarks is smaller than that between quarks in neutrons. Therefore the value of o? in Eq. ffT^ 
for the dimension-2 condensate is expected to be larger in the inner part than in the outer 
part of the neutron star. Because the neutron star is in the form of sphere, the gauge slices 
dehned by Eq. ffT^ have to be spheres and the value o? is expected to be decreasing from 
the center of the neutron star. In order to check out the decreasing pattern of the value of 
dimension-2 condensate, we need to review the change of patterns for various configurations 
of quarks. 

The hrst example to check the variance of dimension-2 condensate is the baryon. If we 
assume that the three quarks are located at the vertices of an equilateral triangle, we can 
draw the equi-(74^) surfaces as in Fig. [3]|^a). In Fig. [3]^b), the variations of values 
along the axes through the center of the triangle are shown as a function of the distance 
from the center. The next example is a tetraquark state which can be a candidate for 
some X particles discovered in the search of (cc) excited states. A typical tetraquark 
conhguration is drawn in Fig. HJ^a) and the variations of dimension-2 condensate are shown 
in Fig. m^b). 

Similarly we can draw the configuration of a deuteron as in Fig. [5]^ a) and the condensate 
values are shown in Fig. [5t^b). The six quarks are taken to be at the corners of a trigonal 
prism and the condensate values are along the 3 lines. The case of ^He can be drawn as in 
Fig. ini In Fig. Et^a), the nine quarks are at the corners of three triangles, and the condensate 
values are shown in Fig. [6]^b). In Fig. [71 another case with the nine quarks on a circle is 
shown. We can find that the condensate values have maximum value at the center except 
for the positions of quarks and decay exponentially just outside of the quark locations. 
The same situation is expected to be valid for the case of quark locations in the shape of 
spherical shell. By combining spherical shells, we can obtain the case of sphere shape quark 
locations. This case is approximately in the same situation as that of a neutron star, but the 
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neutron star contains the region of quark-gluon plasma. For the quarks in the quark-gluon 
plasma region, we cannot apply our nonperturbative formalism in a naive manner but their 
contributions to long range conhning region can be induced from the above simple examples. 
It is expected that in the quark-gluon plasma region the values become saturated due 
to the perturbative interactions but the quark contributions to the {A‘j^) values in the outer 
part of the neutron star could be increased. This situation is shown in Fig. [HI The radius for 
quark-gluon plasma region is represented by Ri and the radius of the whole neutron star is 
designated as Rq. There exist no quarks outside of i?o and because the {Aj^) value extends 
over the space where no baryonic nor nuclear structures exist, it is possible for low energy 
gluons to propagate along the macroscopic spherical surfaces outside of neutron star. These 
gluons are not observed by photons or by neutrinos and therefore can be good candidates 
for dark matter. 

The gluons propagating in free space may have originated from the Big Bang. As the 
universe expands, the energies of gluons decrease below the value of 20 MeV and then they 
cannot be captured in baryons or nuclei. However, they can interact with each other and 
form condensate structures around neutron stars and even out to the scale of a galaxy. If 
the sterile gluons left over after the big-bang nucleosynthesis existed uniformly over the 
universe, they could exist in inter-galactic regions. But somehow if they were gathering 
around galaxies, it is possible for them to move from inter-galactic regions to intra-galactic 
regions. The large scale condensate structures could be formed as nearly spherical shapes 
centered at each centers of mass of galaxies. Since they propagate along the spherical gauge 
slices, the mass effects dehned by Eq. flHOl) could be related to the curvature scalar and we 
get the radius dependences as 


m(r) oc 


(56) 


These dependences are consistent with the ^ variation of the dark matter density in a 


galaxy |38|. For a crude estimation of dark matter density, let’s assume that the average 


energy of gluons propagating in intra-galactic region is approximately 1 MeV fairly below 
the binding energy of deuteron 2.2 MeV. If the average number of gluons per unit volume 
is assumed to be the same as that of Big Bang photons 411 [39|, we have 


10 ^ X 411 ~ 0.4 (GeV/cm^) , 


(57) 


which is just number obtained for the dark matter density in our galaxy. 
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VI. DISCUSSIONS 


In this paper, we reviewed the idea of generalized 9 vacunm and showed various structures 
of gauge hxed slices dehned by the values of dimension-2 condensates. Because the gluons 
have to be quantized on these gauge slices, it is essential to account for the curved structures 
of gauge slices. The ordinary derivatives have to be replaced by covariant ones and then 
from the second order form of the Lagrangian density we get the curvature dependence of 
the gluon propagator. The curvature dependences can be interpreted as mass effects for the 
case of maximum symmetry and these mass effects might be the newly found roles of gluons 
in hadrons and in any strongly interacting systems as dark matter. 

Gluonic contributions to hadronic masses are usually parametrized as constituent quark 
masses which are very large compared with the current quark masses determined by electro- 
magnetically or weakly interacting probes. The differences can be taken as dark matter in 
the same sense of dehning the dark matter of a galaxy by the difference between measured 
and observed masses of a galaxy. As for the large scale conhnement we can consider neutron 
stars as the centers for the large scale structures of gauge slices and the structures can be 
extended even to the whole scale of a galaxy. In general, at the center of a galaxy we expect 
the existence of supermassive black hole and that black hole may serve as the center for the 
galaxial structure of gauge slices. Since the curvature scalar decreases for the outer sphere 
from the center, the mass density of dark matter can be explained systematically from the 
large scale structure of gluon gauge slices. 

Of course, there are many issues to be resolved with more detailed models. First of all, 
the calculation of gluon propagator with more general form of Ricci tensor is necessary. 
This issue is closely related to the interpretation of curvature in terms of the gluonic prop¬ 
erties such as mass or other quantum numbers. In this connection, the argument of gluon 


splitting 


40| could be exploited because we need to separate parallel and transverse spin 


components to follow up the changes along the curved gauge slices. The determination of 
energy spectra of sterile gluons is another issue and the coupling strength of very low energy 
gluons ^ could affect the hnal shape of cosmological dark matter distributions. 

In hadronic scales, the calculations of hadronic masses with appropriate consideration of 
gluonic effects are important issues. Traditionally an appropriate form of conhning poten¬ 
tial is assumed between quarks and the spin-dependent forces 42| are derived by one gluon 
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exchange diagrams in flat space. There exist some divergence problems originating from 
masslessness of the exchanged glnon and we need to resolve the problems with new propa¬ 
gator. From the distribntion of gange slices it is also possible to dednce the form of vacnnm 


distribntion fnnctions which are nsed to define nonpertnrbative qnark propagator 


43|. 
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FIG. 1: (Color online) Gauge fixed slices for the quarks at (—0.2,0.35,0.0), (—0.35,0.15,0.0), 
(-0.25,-0.25,0.0), (-0.2,-0.4,0.0), (0.2,-0.1,0.0), and (0.25,0.05,0.0). 
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FIG. 2: (Color online) Gauge slices for the quarks at (—0.5,0.0,0.0) and (0.5,0.0, 0.0) with (3 = 1.0 
and A: = 1.0. 
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FIG. 3: (Color online) (a) Equi-(^^) surfaces for a baryon with quarks at (0.0, ^,0.0), 
(—0.75,0.0) and (0.75,0.0). (b) values along the axes through the center of the 
triangle. 
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FIG. 4: (Color online) (a) Equi-(74^) surfaces for a tetraquark state with quarks at 
(-^,-0.75,-^), (^,0.0,-^), (-^,0.75,-^) and (0.0,0.0,^). (b) (A^) values along 
the axes through the center of the tetrahedron. 




FIG. 5: (Color online) (a) Equi-(^^) surfaces for a deuteron with quarks at (--^,—0.75,0.75), 
(^, 0.0, 0.75),(-^, 0.75,0.75), (-^,-0.75,-0.75), (^, 0.0,-0.75) and (-^, 0.75,-0.75). (b) 
(A^) values along the x—axis, z—axis and the line with 9 = 90°, (j) = 90°. 


24 













FIG. 6: (Color online) (a) Equi-(^^) surfaces for ^He with quarks at (0.0, ^,0.0), 
(-0.75,-^,0.0), (0.75,-^,0.0), (-0.75, 0.0), (-2.25,-^, 0.0), (-1.50,-^3, 0.0), 

(1.50,—\/3,0.0), (2.25,—^, 0.0) and (0.75, 0.0). (b) (^4^) values along the x—axis, z—axis 

and the line with 6 = 90°, cj) = 324°. 



FIG. 7: (Color online) (a) Equi-(74^) surfaces for a nine quark configuration with 

quarks at (0.0,2.2,0.0), (-1.4,1.7,0.0), (-2.2,0.38,0.0), (-1.9,-1.1,0.0), (-0.75,-2.1,0.0), 
(0.75,-2.1,0.0), (1.9,-1.1,0.0), (2.2,0.38,0.0) and (1.4,1.7,0.0). (b) {A^) values along the two 
axes. 
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FIG. 8: (a) Equi-(yl^) surfaces for a neutron star. The dotted region represents the quark- 

gluon plasma region and the wiggly lines represent the gauge slices outside of the neutron star, 
(b) Expected variation of values from the center of the neutron star. i?j is the radius for 
quark-gluon plasma region and Rq is the radius of the neutron star. 





